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Given an initial quantum state \ipi) and a final quantum state \iPf) in a Hilbert space, there 
exist Hamiltonians H under which \ipj) evolves into \iPf)- Consider the following quantum brachis- 
tochrone problem: Subject to the constraint that the difference between the largest and smallest 
eigenvalues of H is held fixed, which H achieves this transformation in the least time r? For Her- 
mitian Hamiltonians r has a nonzero lower bound. However, among non-Hermitian "PT-symmetric 
Hamiltonians satisfying the same energy constraint, r can be made arbitrarily small without violat- 
ing the time-energy uncertainty principle. This is because for such Hamiltonians the path from \tpi) 
to \iPf) can be made short. The mechanism described here is similar to that in general relativity 
in which the distance between two space-time points can be made small if they are connected by a 
wormhole. This result may have applications in quantum computing. 
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Suppose that one wishes to transform unitarily a state 
in a Hilbert space to a different state by means 
of a Hamiltonian H . In Hermitian quantum mechanics, 
such a transformation requires a nonzero amount of time, 
provided that the difference between the largest and the 
smallest eigenvalues of H is held fixed. However, if we ex- 
tend quantum mechanics into the complex domain while 
keeping the energy eigenvalues real, then under the same 
energy constraint it is possible to achieve such a trans- 
formation in an arbitrarily short time. In this paper we 
demonstrate this by means of simple examples. 

The paper is organized as follows: We first review why 
in Hermitian quantum mechanics there is an unavoidable 
lower bound r on the time required to transform one state 
into another. In particular, we consider the minimum 
time required to flip unitarily a spin-up state into a spin- 
down state. We then summarize briefly how Hermitian 
quantum mechanics can be extended into the complex 
domain while retaining the reality of the energy eigenval- 
ues, the unitarity of time evolution, and the probabilistic 
interpretation. In this complex framework we show how 
a spin-up state can be transformed arbitrarily quickly to 
a spin-down state by a simple non-Hermitian Hamilto- 
nian. Then we discuss the transformation between pairs 
of states by more general complex non-Hermitian Hamil- 
tonians. We make some comments regarding possible 
experimental consequences of these ideas. 

In Hermitian quantum mechanics how does one achieve 
the transformation \ipj) — ► \?Pf) = e~ lHt l h \ipi) in the 
shortest time t = r? Since r is the minimum of all possi- 
ble evolution times t, the Hamiltonian associated with r 
is the "quantum brachistochrone" Finding the opti- 
mal evolution time requires only the solution to a much 
simpler problem, namely, finding the optimal evolution 
time for the 2x2 matrix Hamiltonians acting in the two- 
dimensional subspace spanned by and \i/jf) EJ- 



To solve the Hermitian version of the two-dimensional 
quantum brachistochrone problem one can choose the ba- 
sis so that the initial and final states are given by 
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with \a\ 2 + \b\ 2 = 1. The most general 2x2 Hermitian 
Hamiltonian has the form 
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where the four parameters r, 



eigenvalue constraint 



u, and 8 are real. 
= uj reads 



u? = (s - u) 2 + 4r 2 . 
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The Hamiltonian H in J3J can be expressed in terms of 
the Pauli matrices as H = i(s + u)l + ^ujcr-n, where 
n = h ( rcos ^j' rsm ^j ^~T~) i s a umt vector. Using the 
matrix identity exp(i<^<r-n) = cos^l + i shx cf> <r ■ n, the 
relation \i/jf) = e~ lHr / h \-ipi) takes the form 



a 



' cos (i) - 1 

UJ 
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From the second component of J3J we obtain \b\ = 
psin(|^), which gives the time required to transform 



the initial state 
relation over all r 



t = 2j arcsin (^r). We optimize this 
> 0, keeping in mind that © gives a 



maximum value of \lo for r, achieved when s = u. The 
optimal time is thus 



t = — arcsm \ o\. 



(5) 
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Note that if a = and b = 1 we have r = 2irh/oj for the 
smallest time required to transform ( * ) to the orthogonal 
state ( " ) . This value of r is called the passage time Q ■ 
For general a and b, at the optimal time r we have 



a = e -W?yi _ |ft]2 and 6 



|6|e* & , which sat- 



isfies the condition |a| 2 + \b\ 2 = 1 that the norm of the 
state does not change under unitary time evolution. The 
parameters s and 9 are determined by the phases of a 
and b. Writing a = |a|e 4arg ( Q ) and b = \b\e l ^ b \ we find 
that the optimal Hamiltonian is 



warg(a) 



w -i[arg(b)-arg(a)-§] 
4 e 



rr j 2 arcsin |6 

~~ I w t[arg(b)-arg(o)-f] gaxg(g) 
\ 4 e ~ 2 arcsin 1 6 
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Since the overall phase of |^>f) is not physically relevant, 
the quantity arg(o), for example, may be chosen arbitrar- 
ily and without loss of generality we may assume that it 
is 0. We are free to choose arg(a) because there is no 
absolute energy in quantum mechanics; one can add a 
constant to the eigenvalues of the Hamiltonian without 
altering the physics. Equivalently, this means that the 
value of t cannot depend on the trace s + u of H. 

Interpreting the result for r in JSJ requires care be- 
cause while this equation resembles the time-energy 
uncertainty principle, it is merely the statement that 
ratextime= distance. The constraint J3J on H is equiv- 
alent to placing a bound on the standard deviation 
AH of the Hamiltonian, where AH in a normalized 
state \ip) is given by (AH) 2 = (^\H 2 \^) - (i/j\H\ij} 2 . 
The maximum value of AH is u>/2. According to the 
Anandan-Aharonov relation p|, the speed of evolution 
of a quantum state is given by AH. The distance be- 
tween the initial state |^>j) and the final state \4>f) is 
5 = 2arccos(|( , 0_F|'!/'/)|)- Thus, the shortest time r to 



achieve the evolution from \ipi) to \tJjf) = e 



-iHr/n 



\ipi) is 



bounded below because the speed is bounded above while 
the distance is held fixed. The Hamiltonian H that re- 
alizes the shortest time evolution can be understood as 
follows: The standard deviation AH of the Hamiltonian 
in @ is r. Since AH is bounded by w/2, to maximize the 
speed of evolution (and minimize the time of evolution) 
we choose r — uo/2. 

The objective of this paper is to perform the same 
optimization for complex non-Hcrmitian Hamiltonians 
having PT symmetry. There are infinitely many PT- 
symmetric complex non-Hermitian Hamiltonians whose 
eigenvalues are real and bounded below. Here, P is the 
parity reflection operator and T is the time reversal op- 
erator. A one-parameter family of such Hamiltonians 
that has been investigated intensively 0, is given by 
H = p 2 + x 2 (ix) e , where e > 0. Although this Hamilto- 
nian is not Hcrmitian in the usual Dirac sense, where Her- 
mitian adjoint consists of complex conjugation and ma- 
trix transposition, H defines a unitary theory of quantum 
mechanics 00- This is because H is self-adjoint with 
respect to a new inner product that is different from the 



Dirac inner product of conventional quantum mechanics. 

This new inner product is expressed in terms of a linear 
operator C that satisfies three equations Q : 



C 2 = 1, [C, H] = 0, and [C, PT] = 0. 



(7) 



For any given H we can, in principle, calculate C by solv- 
ing the three equations in (JZj). We then define an inner 
product in terms of CPT conjugation. Thus, in a PT- 
symmetric quantum theory the inner product is specified 
dynamically in terms of the Hamiltonian. Furthermore, 
the time-evolution operator e ~ lHt / h j s unitary (norm pre- 
serving) because H commutes with CVT . 

It has been shown that for any 'PT-symmetric Hamil- 
tonian having real eigenvalues there exists an equivalent 
Hermitian Hamiltonian [T5| ■ The argument goes as fol- 
lows: The Hermitian operator CV is positive and can 
thus be written as CV = e®. By means of the similarity 
transformation H = e~^/ 2 He^/ 2 one can construct the 
corresponding Hcrmitian Hamiltonian H, which is equiv- 
alent to H in the sense that it has the same eigenvalues. 
The states in the PT-symmetric theory are mapped by 
the operator e~®/ 2 to corresponding states in the Her- 
mitian theory. But, since this operator does not keep 
the states in the same Hilbert space, relative properties 
of states can be changed. For example, the overlap dis- 
tance between two states does not remain constant in 
the original Hilbert space. In this paper we exploit this 
property of the e~^/ 2 transformation to circumvent the 
Hermitian limit on t. 

We now show how to solve the 'PT-symmetric brachis- 
tochrone problem for a simple non-Hermitian two- 
dimensional matrix Hamiltonian of the form 



H = 



re 
•s 



re 
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(This Hamiltonian was examined in detail in Ref. jjj •) 
To show that H in JHJ is PT symmetric, we let T be the 
operation of complex conjugation and V be given by 



V = 



1 

1 



(9) 



The eigenvalues E± 



± Vs 



r 2 sin 2 9 of H in 



(JSJ) are real provided that s 2 > r 2 sin 2 9. This inequality 
defines the region of unbroken PT symmetry. 
The unnormalized eigenstates of H are 



\E + ) = 



e ia/2 
e — ia/2 



\E_) = 



ie~ la l 2 
-ie ia l 2 



(10) 



where the real parameter a is defined by sin a — 
(r/s)sui9. The operator C satisfying the conditions in 
Q is given by 



C = 



1 



cos a 



1 

-i sin a 



(11) 
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By using @ and we calculate that the CPT norms 
of both eigenstates in (jTU)l are \/2 cos a. 

Following the procedure used for Hermitian Hamilto- 
nians, we rewrite H in JSJ in the form H = (r cos 9)1 + 
iwer-n, where n = ^(s, 0, irsm9) is a unit vector and 
the squared difference between the energy eigenvalues is 



u? = 4s 2 - At 2 sin 2 ( 



(12) 



The positivity of uj 2 is ensured by the condition of un- 
broken PT symmetry. This equation emphasizes the key 
difference between Hermitian and PT-symmetric Hamil- 
tonians: The corresponding equation J2JI for a Hermi- 
tian matrix Hamiltonian has a sum of squares, while this 
equation for uj 2 has a difference of squares. Thus, Hermi- 
tian Hamiltonians exhibit elliptic behavior, which leads 
to a nonzero lower bound for the optimal time r. The hy- 
perbolic nature of (fL2"|> allows r to approach zero because, 
as we will see, the matrix elements of the PT-symmetric 
Hamiltonian can be made large without violating the en- 
ergy constraint E + — E- = uj. 

The 'PT-symmetric analog of the evolution equation 
(@J) is given by 



-iHt/h 



-itr cos 9 j 



COS a 



n ( co s (i " 



-i sm 



(13) 



We apply this result to the pair of vectors examined in the 
Hermitian case: = (q) and \i^f) = (])• (Note that 
these vectors are not orthogonal with respect to the CPT 
inner product.) From Ijl3|l we see that the evolution time 
to reach \"4>f) from is t = (2a + it)h/uj. Optimizing 
this result over allowable values for a, we see that as a 
approaches —\k the optimal time r tends to zero. 

Note that in the limit a — > — i-7r we get cos a — > 0. 
However, in terms of the variable a the energy constraint 
(I12J1 becomes ui 2 = 4s 2 cos 2 a. Since uj is held fixed, in or- 
der to have a approach — ^7r we must require that s>l. 
It then follows from the relation sin a = (r/s) sin 9 that 
\r\ ~ |s|, so we must also require that r 1. Evidently, 
in order to make r <C 1 the matrix elements of the PT- 
symmetric Hamiltonian (JHJ must be large. 

The result demonstrated here does not violate the un- 
certainty principle. Indeed, a Hermitian Hamiltonian 
and a non-Hermitian PT-symmetric Hamiltonian both 
share the properties that (i) the passage time is given by 
2nH/cu, and (ii) AH < ui/2. The key difference is that a 
pair of states such as ( J ) and ( 1 ) are orthogonal in a Her- 
mitian theory but have the separation 5 = ir — 2\a\ in the 
PT-symmetric theory. This is because the Hilbert space 
metric of the PT-symmetric quantum theory depends on 
the Hamiltonian. As a consequence, it is possible to set 
the parameter a to create a wormhole-like effect in the 
Hilbert space 

A gedanken experiment to realize this effect in a lab- 
oratory might work as follows: We use a Stern-Gerlach 



filter to create a beam of spin-up electrons. The beam 
then passes through a 'black box' containing a device 
governed by a PT-symmetric Hamiltonian that flips the 
spins unitarily in a very short time. The outgoing beam 
then enters a second Stern-Gerlach device that verifies 
that the electrons are now in spin-down states. In ef- 
fect, the black-box device is applying a magnetic field in 
the complex direction (s, 0, ir sin 9). If the field strength 
is sufficiently strong, then spins can be flipped unitar- 
ily in virtually no time because the complex path join- 
ing these two states is arbitrary short without violating 
the energy constraint. The arbitrarily short alternative 
complex pathway from an up state to a down state, as 
illustrated by this thought experiment, is reminiscent of 
the short alternative distance between two widely sepa- 
rated space-time points as measured through a wormholc 
in general relativity |l2( • 

The PT-symmetric Hamiltonian (jSJ) used in the fore- 
going illustrative example contains only three arbitrary 
real parameters, which are not sufficient to allow the ini- 
tial state \ipi) = (J) to evolve into any final state \iPf}- 
Indeed, it follows from (|13fl that (i) can only evolve into 
(^) if the relative phase of a and b is ±5^. Therefore, 
we introduce the more general four-real-parameter PT- 
symmetric Hamiltonian 



H 



x + (z + iy) 

z 

tan 7 



7 iy tan 7 

tan 7 « ' 

?j/tan7 x — (z + iy) 



(14) 



which is associated with a more general definition of par- 
ity reflection P than that used in © : 



P 



sm 7 cos 7 
cos 7 — sin 7 



(15) 



We retain the same definition for T, namely, that T per- 
forms complex conjugation. 

As before, we express H in the form H — xl + ^cr n, 
where in this case n = ^(^/tan7 — iyta,nj,0,z + iy). 
The operator C — cr ■ n is given by 



C = 



z + iy 



tan 7 



uj \ iutan7 

v tan 7 y ' 



iy tan 7 
-z - iy 



(16) 



which, along with l|15f) . allows us to define the inner prod- 
uct with respect to which the Hamiltonian (f 141) becomes 
self- adjoint. The energy constraint E+ — E- = uj again 
takes a hyperbolic form: 



uj = Az esc 7 — 4y sec 7. 



(17) 



For this more general Hamiltonian, the initial state (J) 
evolves as follows: 



-iHt/h 



1 



cos(i) + ^sin(i)^sin(g) 



, 2h, 



2y 



tan 7 sin g 



2.3 



uj tan 7 



.(18) 
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The time evolution preserves the CPT norm of the ini- 

/ 9 \ 1/2 

tial state, which is ( . . We therefore choose the 

1 y u sin 7 I 

general form of the final state to be 



IV-'j- 



ue 

u) sin 7 \ 



(19) 



where u, v, A, and £ are real parameters. 

We now introduce dimensionless variables X — 2x/ui, 
Y = 2y/uj, Z = 2z/ojsm^f, and T = uit/2h, as well as the 
shifted variable B = A + xt/h. By identifying the right 
sides of (I18fl and i|19|) . we obtain 



cos T + Y sin T = y/Zu cos B, 
ytan7sinT= -y/Zv cos(B + £), 
\/l?sin7sinT = — usinB, 
"\/Zcos7sinT = — usin(_B + £). 



(20) 



The energy constraint in l|17ll takes the hyperbolic form 

1 = Z 2 - Y 2 sec 2 7. (21) 

The condition that the norm of the initial vector be 
preserved under time evolution imposes the requirement 
that u 2 + v 2 + 2itvFsin£/Zcos7 = 1, which can be de- 
rived from the five equations in (|20|l and l|21|) . 

The generic problem is now to pick a final vector; that 
is, to choose the parameters A, u, v, and £ in 119|) . and 
then to solve (|2C)|) and (|21|l to determine the parameters 
X, Y, Z, and 7 for which \tpj) reaches \iPf) under the 
Hamiltonian in l|14|) . We then must find the smallest 
value of T for which the transformation is realized. 

To illustrate the procedure, we consider the example 
u = v and solve the five simultaneous equations in the 
form of Laurent series valid for small T. The result is: 



x = (h-U~ A ) t^ + oco, 

Y = -uT- 1 cos(±7r - ±£) + O(T), 
Z = -uT- 1 + 0(T), 
7 = i7r-i£ + 0(T 2 ). 



(22) 



Note that the parameter T may be taken arbitrarily small 
and thus the initial vector evolves into the final vector in 
an optimal time that is arbitrarily small. Of course, the 
matrix elements of the Hamiltonian become large in this 
limit, but this is possible because the energy constraint 
in l|17l) is hyperbolic in character. 

We conclude by remarking that the results established 
here provide the possibility of performing experiments 
that definitively distinguish between Hermitian and VT- 
symmetric Hamiltonians. If practical implementation of 
complex 'PT-symmetric Hamiltonians were proved feasi- 
ble, then the identification of the optimal unitary trans- 
formation would be particularly important in the design 
and implementation of fast quantum communication and 
computation algorithms (cf. ^|). 



Of course, the wormhole-like effect we have discussed 
here can only be realized if we can switch rapidly between 
Hermitian and PT-symmetric Hamiltonians by means of 
similarity transformations. It is conceivable that there 
is a sort of quantum protection mechanism that places 
a lower bound on the time required to switch Hilbcrt 
spaces. If so, this would limit the applicability of a 
Hilbert-space wormhole to improve quantum algorithms. 
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